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Abstract 



1 Introduction 

Recent advances in nanofabrication technology have made possible to cre- 
ate low dimensional structures [ 1-4 and references therein] such as quantum 
wells, wires and dots. In almost most of problems concerning such struc- 
tures the eigenspectrum of the 2D- Schrodinger equation is investigated. For 
example, the 2D hydrogenic energy levels in a constant magnetic field of ar- 
bitrary strength has been a subject of numerous theoretical and experimental 
investigations [1-14]. 

Taut [6] has found analytical solutions for this problem. However, no solu- 
tions for nodeless states were found. Martin et al [10] have used a two - point 
quasifr actional approximation and found excellent interpolation between the 
weak and strong magnetic field perturbation expansions. MacDonald and 
Ritchie [8] have however used a two - point Fade approximation. The curves 
obtained from different Fade approximates are very different. No regular 
pattern appears and the results become unreliable. More detailed discussion 
can be found in ref [10]. Zhu et al [9] have used a power series expansion 
method. Their results compared excellently with the direct numerical inte- 
gration [12] and perturbation [8] methods in the weak field regime. However, 
neither Zhu's nor direct numerical integration results seem to approach the 
perturbation results at the strong field limit. Villalba and Fino [11] have 
used a finite - difference scheme with a linear mesh of up to 2000 points and 
failed to provide a good estimation of the ground state for the 2D hydrogen 
atom. On the other hand, their variational solution obtained using hydrogen 
basis is in good agreement with that of Martin et al [10] at the weak field 
limit. Using the oscillator basis, their results were in good agreement at the 



strong field limit. Their solutions could not therefore provide information 
about how the energy shifts in the intermediate range of the magnetic field. 
Mustafa [5] and Quiroga et al [7] have used the shifted N - expansion tech- 
nique (SLNT) and their results were in good agreements with those of Martin 
et al [10]. 

In this work, section two considers the 2D Shcrodinger equation with 
an arbitrary cylindrically symmetric potential. By introducing a technique, 
which we will call Pseudoperturbative Shifted - / Expansion Technique (PSLET), 
we obtain analytical expressions for both eigenenergies and eigenfunctions for 
the 2D Schrodinger equation. In section 3, we test our method for both 2D 
Coulomb and harmonic oscillator potentials, and the results found to be ex- 
act. We also consider the 2D hydrogenic energy levels in a constant magnetic 
field of arbitrary strength. PSLET results are shown to compare excellently 
with the perturbation [8], direct numerical integration [12], and series expan- 
sion [9] methods. They are also shown to be in exact agreements with the 
numerical solution of the 2D Schrodinger equation [10]. We conclude with 
remarkes in section 4. 

2 Theory of 2D-PSLET 

The Schrodinger equation for an arbitrary cylindrically symmetric potential 
V{p) { in h = 2m = 1 units) is 



d' 4/2 - 1 ^^^ ■ 
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^ip)=E^ip), (1) 



where p'^ = x'^ + y'^ , I = \m\ and m is the magnetic quantum number. If / is 
shifted through the relation I = I — P, Eq.(l) becomes 
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Where Q is a scale to be determined below and set equal to /^ at the end of 
the calculations, and f3 is a suitable shift to be determined and is introduced 
partly to avoid the trivial case / = 0. 

The systematic procedure of PSLET starts with making use of Taylor's 
theorem and expanding (2) about an arbitrary point ( for now) on the p - 
axis. It is convenient then to transform the coordinates in (2) via the relation 



x = l^/^{p-po)/po, 
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where Po is our arbitrary point. Expansions about a; = yield 
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where 
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It is convenient to expand E as 



E= ^ E(")/-". 



(7) 



n=-2 



where the prime of V{po) denotes derivative with respect to po- Equation (6) 
when compared to Schrodinger equation for one - dimensional anharmonic 
oscillator 
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where P{x) is a perturbation term and So is a constant, implies 



£o = I 
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and 
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Equations (10) and (11) yield 
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Here po is chosen to minimize E^ ^\ i. e. 



= and — — :; — > 0, 



dpo 
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which in turn gives, with / = ^/Q 
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The shifting parameter j3 is determined by choosing the next leading cor- 
rection to the energy eigenvalue, E^~'^\ to vanish. This choice is physically 
motivated by requiring the agreement between 2D-PSLET eigenvalues and 
eigenfunctions and the exact analytical eigenvalues and eigenf unctions for 
both the Coulomb and harmonic oscillator potentials. Hence 
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where 
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Then equation (6) becomes 
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where 
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Equation (5) thus becomes 
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For nodeless unnorinalized wave functions Up = and 



^o(a;(p)) = ea;p(f/(a;)), 



(26) 



which when substituted in equation (25) yields 
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where prime of U{x) denotes derivative with respect to x. It is evident that 
this equation admits solution of the form 
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Substituting equations (28) into equation (27) implies 
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where primes of U^'^\x) and G''^"'^(a;) denote derivatives with respect to x. 
Equating terms of same order in / one obtains 



_[;7(o)'+ ^(0)^(0)] ^^(0)^0^ 



(32) 



f/(°) (a;) = A,o ; ^1,0 =-uj/2 

Integration over dx yields 

U^^\x) = -ujx/2 

Similarly, 

_ [f/(i)' + G-W'] - 2f/(o)f/(i) - 2f/(o)G'(°) + t;(i) = 0, 

f/«(a;) = 0, 
9 



(33) 



(34) 



(35) 
(36) 



G'(°)(^) = Co,o + Ci,oa;2, (37) 



Ci,o = -— , (38) 

w 



Co,o = ^{2Ci,o-w), (39) 



^--^ + |^^(^°)' (^°) 
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f/(2)(a;) = Dl,2a: + /^2,2a:^ (42) 

G^^\x) = 0, (43) 

^2,2 = -(Ci2o-52) (44) 

w 



Z?i,2 = ^(3/^2,2 + 2Co,oCi,o - 6/5), (45) 



^-^+2!^^^^°^' ^''^ 



10 



A(°) = -(A,2 + 0. (47) 



... and so on. Clearly one can calculate the energy eigenvalues and the 
eigenfunctions from the knowledge of Cn,m and D„ „ in a hierarchical man- 
ner. However, it is for the convenience of this study to conclude the procedure 
here and give the energy eigenvalues and eigenfunctions as 

E = E^-^^P + E(o) + E^'H-' + E^-^ + 0{l~^), (48) 



u\x{p)) = f/(o) + [G(o) + f/«]ri/2 + [G« + f/(2)]ri 

+ [G'(5)+f/(6)]r3 + 0(r^/2). (49) 



From which we find 



'^q{x{p)) = expij U\x)dx\ . (50) 



3 Applications, results and discussions 

In this section the above analytical expressions of 2D - PSLET are investi- 
gated through the 2D Coulomb, harmonic oscillator and the hybrid of the 
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two potentials. 

3.1 The Coulomb potential V{p) = -2/p 
Following 2D - PSLET theory one obtains, for nodeless states, 

E^-^^? = -l ■ l = \m\ + l/2 , p, = /^, (51) 

^(0) ^ ^{1)1-1 ^ ^(2)J-2 = . . . = 0, (52) 



ui.) = -iy+i(y-y- + y--y-+y--y- + y--y-], (53) 

^ ' ^ V 2 3 4 5 6 7 8/' ^ ' 



where y = xl^^^"^. It is obvious that the second term in equation (85) is the 
infinite geometric series expansion for /n(l +y)'. Equation (85) can thus be 
approximated by 

U{x) ~ -ly + Zn(l + y)\ (54) 

which in turn implies that 

^o,™(p) = [-) e^e-'^'P°. (55) 

Hence the nodeless radial parts of the wave functions are 

RoAp)=Np'-^e''p'P% (56) 
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where N is the normahzation constant given by the relation 

7V = M^^ (57) 



(2()! 

Equation (88) evidently gives the exact expressions for the normalized node- 
less radial parts of the wave functions. 

Finally, the exact energy eigenvalues for any nodeless orbital state are 

Eo,m = -{\m\ + l/2)-2. (58) 



3.2 The harmonic oscillator potential V{p) = 7^/9^/4 
For this potential, 2D - PSLET procedure yields, for nodeless states, 

E^-^^F = -fl ; l = \m\ + l, (59) 

^(0) ^ ^{1)1-1 ^ e{2)J-2 = . . . = 0, (60) 



../ X 1 / y\ y^ y* , y^ y^\ 

U[x) = — \y 

^ ' 2V2 3 4 5 6/ 

,j[ y\ y^ y^ , y^ y\ y'^ y^\ 
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Obviously the terms in brackets are the infinite geometric series expansions 
for ln{l + y). Equation (93) can thus be approximated by 

U{x) - ln{l + i/)-^/2 ^ i^(^i ^ yY _iy_ i!L^ (62) 



which in turn imphes that 



1-1/2 



^oAp)={y] e-'^'/' ■ Po = plh. (63) 



Hence, the nodeless radial parts of the wave functions are 

R,Ap) = Nt'/'e~^P''\ (64) 

where the normalization constant N is given through the relation 



N' = -^ '^" ' -J^. (65) 

1-3 -5 (2/-2)V27r ^ ^ 



Equation (96) clearly yields the exact expressions for the normalized node- 
less radial parts of the wave functions. The exact energy eigenvalues for any 
nodeless state are given by 

Eo,m = l{\m\ + 1). (66) 
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3.3 A hybrid of Coulomb and oscillator potentials 

Perhaps the most interesting form of such a hybrid model is 

V{p) = m7 - - + ^, (67) 

where m is the magnetic quantum number. This potential (99) describes a 
2D electron gas in the x-y plane in the presence of a hydrogenic potential, 
representing the interaction between a conduction electron and a donor im- 
purity, in a magnetic field in the z - direction. It also resembles the potential 
model that describes a donor impurity in a single quantum dot moving in a 
magnetic field in the z - direction [4-14]. 

It is obvious from equation (17), along with (18) and (19), that one can 
hardly find an analytical solution for po-, in terms of 7 and m, using the 
potential in (99). One has therefore to appeal to numerical techniques to 
solve for p^ for each 7 and m. Once po is found, the energy eigenvalues 
as well as eigenfunctions can be obtained through the 2D-PSLET theory 
described in section 2. 

In tables 1-4 we list the energy eigenvalues in such a way that the con- 
tribution of each energy correction is made clear. The tables show that 
2D-PSLET results are rapidly convergent. It is also evident that the energy 
eigenvalues for states with larger \in\ converge more rapidly than those with 
smaller \m\. Such a tendency was obvious from the very moment of the in- 
vention of the deceptive perturbation parameter 1//, / = \m\ — (3 and (3 is 
always negative. 

It could be interesting to know that the computation time for all po's 
needed for the entries in tables 1-4 is less than 30 sec, and for each entry is 
at most 20 sec, including the eigenfunctions for each of them, po is computed 
using EUREKA and the eigenvalues as well as eigenfunctions, in terms of p, 
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are computed using REDUCE 3.4 on a standard Pentium PC. 

In figure 1 ( Available from authour upon request ) , 2D-PSLET results 
(solid curve) for the IS-state compare excellently with the others in the weak 
magnetic field regime. In the strong field regime, 2D-PSLET results fall in 
between the perturbation (small dashes) [8] and direct numerical integration 
(solid curve connecting solid circles) [12], on which the series expansion re- 
sults (solid circles on the solid curve) [9] are located. SLET and SLNT [3] 
results are presented by the long dashed curve. 2D-PSLET, SLET and SLNT 
results are, however, unique in their tendency to approach the strong and the 
weak magnetic field perturbation results [8]. The perturbation coupling con- 
stants were appropriately defined in these regimes. Likewise, we believe, it 
should be the tendency of the results of any approximation technique. 

Figures 2-4 ( Available from authour upon request ) show that 2D-PSLET 
results (empty squares) are in exact agreements with those of the numerical 
solutions of the Schrodinger equation [10]. The best fit line of 2D-PSLET 
results exactly overlaps with that of the lowest two - point quasifr actional 
approximation [10]. 

It should be mentioned that, for all states considered above, 2D-PSLET 
results for the energies, excluding the third - order correction, are exactly the 
same as those of SLET [3] and SLNT [3,5]. It is only because of the third - 
order correction, E^'^^P, that the 2D-PSLET results compared better, with 
the numerical results, than those of SLET and SLNT. 

4 Conclusions and Remarks 

In this work, the psuedoperturbative shifted - / expansion technique (PSLET) 
was introduced to find nodeless states of 2D Schrodinger equation with ar- 
bitrary cylindrically symmetric potentials. Exact energy eigenvalues and 
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eigenfunctions for 2D Coulomb and harmonic oscillator potentials were re- 
produced. Also, exact energy eigenvalues, compared to numerical ones [10], 
were obtained for the hybrid of the 2D Coulomb and oscillator potentials. 
The accuracy and rapid convergence of 2D-PSLET results are satisfactory 
and fascinating. The analytical results, tables, and figures clearly bear this 
out. 

Finally, some observations concerning the attendant 2D-PSLET are in or- 
der. It is highly accurate and rapidly convergent, thus efficient with respect 
to computation time. Within the same procedure, it produces both eigenval- 
ues and eigenfunctions. It puts no constraints on the coupling constants of 
the potential involved. It is to be understood as being an expansion through 
any existing quantum number in the centrifugal - like term of any Scrodinger 
- like equation, equation (1). 
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Table 1: PSLET results for IS-state energies (in effective Rydberg units) 
of the 2D hydrogenic levels in a magnetic field , where ENq = PE^~'^\ 
ENi = ENo + E(°), EN2 = ENi + E(^)/r, and EN^ = EN2 + E'^^^/P. 



1 


EiVo 


ENi 


EN2 


EN^ 





-4.0 


-4.0 


-4.0 


-4.0 


I 


-3.936821 


-3.907949 


-3.910053 


-3.910538 


2 


-3.738814 


-3.651631 


-3.677083 


-3.673240 


3 


-3.378562 


-3.284620 


-3.356622 


-3.314095 


4 


-2.810985 


-2.864030 


-2.934645 


-2.900042 


5 


-2.024429 


-2.380567 


-2.401652 


-2.496601 


6 


-1.141379 


-1.784237 


-1.840713 


-1.971307 


7 


-0.286914 


-1.115087 


-1.264520 


-1.373630 


8 


0.525016 


-0.421911 


-0.660921 


-0.761824 


9 


1.311535 


0.279818 


-0.032861 


-0.139894 


fO 


2.085744 


0.987287 


0.614648 


0.494518 


f2 


3.623940 


2.419691 


1.954631 


1.802345 


20 


9.888892 


8.387814 


7.701685 


7.438987 


28 


16.38730 


14.65978 


13.83081 


13.48993 


36 


23.06128 


21.13973 


20.19640 


19.79357 


40 


26.44867 


24.43836 


23.44415 


23.01409 
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Table 2: PSLET results for IS-state energies (in effective Rydberg units) 
of the 2D hydrogenic levels in a magnetic field , where ENq = PE^~'^\ 
ENi = ENo + E(°), EN2 = ENi + E'^'^^/l EN3 = EN2 + E^^^/?, and 7' 
= 7/(1 + 7). 



7 


ENo 


ENi 


EN2 


EN^ 


0.0 


-4.0 


-4.0 


-4.0 


-4.0 


0.1 


-3.999228 


-3.998843 


-3.998843 


-3.998830 


0.2 


-3.996091 


-3.994148 


-3.994156 


-3.994162 


0.3 


-3.988498 


-3.982838 


-3.982914 


-3.982940 


0.4 


-3.972089 


-3.958670 


-3.959106 


-3.959233 


0.5 


-3.936821 


-3.907949 


-3.910053 


-3.910542 


0.6 


-3.855896 


-3.797597 


-3.807128 


-3.807128 


0.7 


-3.638726 


-3.538277 


-3.578849 


-3.566445 


0.8 


-2.810985 


-2.864030 


-2.934645 


-2.900042 
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Table 3: PSLET results for 2P - state energies (in effective Rydberg units) 
of the 2D hydrogenic levels in a magnetic field , where ENq = PE^~'^\ 
ENi = ENo + E^°\ EN2 = ENi + E^'^^l, EN3 = EN2 + E^^^P, and 7' 
= 7/(1 + 7). 



7 


ENo 


ENi 


EN2 


EN, 


0.0 


-4/9 


-4/9 


-4/9 


-4/9 


0.1 


-0.525152 


-0.523702 


-0.523991 


-0.523943 


0.2 


-0.556955 


-0.562580 


-0.562780 


-0.562971 


0.3 


-0.542085 


-0.561985 


-0.563083 


-0.563600 


0.4 


-0.481468 


-0.514474 


-0.518498 


-0.519116 


0.5 


-0.355319 


-0.400522 


-0.408080 


-0.409164 


0.6 


-0.113944 


-0.172226 


-0.183710 


-0.185583 


0.7 


0.366646 


0.292098 


0.275860 


0.272867 


0.8 


1.478312 


1.379563 


1.356515 


1.351825 
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Table 4: PSLET results for 3D - state energies (in effective Rydberg units) 
of the 2D hydrogenic levels in a magnetic field , where ENq = PE^~'^\ 
ENi = ENo + EW, EN2 = ENi + E(i)/r, EN3 = EN2 + E^^yP, and 7' 
= 7/(1 + 7). 



7 


ENo 


ENi 


EN2 


EN, 


0.0 


-4/25 


-4/25 


-A/25 


-A/25 


0.1 


-0.257041 


-0.260346 


-0.260476 


-0.260508 


0.2 


-0.264976 


-0.272340 


-0.273106 


-0.273176 


0.3 


-0.225665 


-0.236222 


-0.237588 


-0.237735 


0.4 


-0.134531 


-0.148204 


-0.150156 


-0.150392 


0.5 


0.031455 


0.014366 


0.011785 


0.011447 


0.6 


0.325970 


0.304776 


0.301450 


0.300990 


0.7 


0.881999 


0.855342 


0.851037 


0.850414 


0.8 


2.115430 


2.080314 


2.074507 


2.073635 
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